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1) Label root voxel with integer 0; 

2) Construct a queue and line up the root in the queue; 

3) If( There is at least one voxel in the queue ) 

Serve the voxel x on the top of the queue: 
For( each of x's 26-connected neighbor voxel y) { 
If( y in the volume and has not been labeled yet ) { 
Line up the y in the queue; 

/* label the voxel y */ 
Set dist = 999999 ; 

For( each of y's 26-connected neighbor voxel z ) { 

If( z in the volume and has been already labeled with an integer of n z ) { 

d z =n z +d(y,z)\ 

where d(y,z) is 10, 14 or 17 if the Euclidean distance between y and z is 
, respectively; 

} 

If( dist > d z ) { 

label y with integer dist ; 
dist = d z ; 

} 

} 

} 

} 

x leaves the queue; 

} 

Else { 

end of calculating the distance map. 

} 



